TAUP 2307-95 
CBPF NF-085/95 
UCI 95-34 
December 1995 



A Screened BFKL Interpretation of F<i in 
the Exceedingly Small x Limit 

E. GOTSMAN a ' 6 '5 E. M. LEVIN^'Q and U. MA0R a,c 'Q 

a School of Physics and Astronomy 
Raymond and Beverly Sackler Faculty of Exact Science 
Tel Aviv University, Tel Aviv, 69978, Israel 

b Physics Department, University of California 
Irvine, CA 92717, USA 

c LAFEX, Centro Brasileiro de Pesquisas Ffsicas (CNPq) 
Rua Dr. Xavier Sigaud 150, 22290 - 180, Rio de Janiero, RJ, Brasil 

d Theory Department, St. Petersburg Nuclear Physics Institute 
188350, St. Petersburg, Gatchina, Russia 



Abstract 

In this letter we show that the behaviour of F2, at very small xb, 
agrees with the behaviour expected from the BFKL evolution equa- 
tion, when screening corrections are included. We obtain a description 
which is consistent with the data, however, we require the screening 
corrections to be relatively large (about a quarter of the total DIS 
cross section). The relation between the screening corrections and the 
diffr active DIS cross section is discussed. 
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In this letter we discuss the dependence of F 2 , the proton structure func- 
tion, on W, the cm. energy, at very small xb- We are motivated by the 
recently published data taken at HERA by the ZEUS [|TJ and HI || collab- 
orations which are shown in Fig.l. Our goal is to extract new information 
from the experimental data on the deep inelastic scattering (DIS) process, 
in the region of very small Xb- We will show that the BFKL evolution equa- 
tion (the BFKL Pomeron) ||, including screening (shadowing) corrections @], 
provides a good reproduction of the observed data. 

We list first the main qualitative properties of the behaviour of F 2 (W, Q 2 ), 
as observed by the two experimental groups at HERA and shown in Fig.l. 

1) For W values below 130-150 GeV, the measured data points cluster in a 
narrow linear band. i.e. F 2 (W, Q 2 ) is approximately linear in W and has a 
weak dependence on Q 2 . 

2) For higher values of W, F 2 {W, Q 2 ) is dependent on Q 2 . The high Q 2 data 
differs from the lower Q 2 data which seems to reach a local plateau, resulting 
in a F 2 which is almost constant as a function of W. 



The features of the data at low W are compatible with a dominance of 
the BFKL Pomeron|| in the small xb domain of F 2 (xb, Q 2 )- This is readily 
seen when we write || the BFKL generated structure function 
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where Go denotes the unknown normalization of the gluon distribution at 
Xb ~ 1- The value of ujq is given by[[| luq = 4 In 2 .In the following 

we assume that t^o = 0.5, which corresponds to a resonable value of a s . At 
small values of xb, W 2 = — , so we can rewrite Eq.(l) in the form 
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The above expression reproduces the qualitative F 2 features listed for the 
lower W, but fails to reproduce the required high energy characteristics. To 
improve the "theoretical" behaviour of F 2 at higher values of W, we introduce 
the shadowing correction illustrated in Fig. 2. The DIS structure function 
can be represented as 

F 2 (x B ,Q 2 ) = F 2 BFKL (x B ,Q 2 ) + AF 2 (x B ,Q 2 ) (3) 

where AF 2 represents the changes in the BFKL structure function which 
result from screening. 

We shall elaborate on the details of the screened diagram calculation, 
shown in Fig. 2, later. Our discussion is based on the main results of Ref.|| 
which are summarized as follows: The dominant contribution of interest 
comes from the exchange of two BFKL ladders (Pomerons), while the upper 
blob of Fig.2 is suitably given by the GLAP|[j DIS structure function. The 
integration over xp (see notations in Fig.2) results in the contribution 
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where we have absorbed all nonperturbative QCD contributions in the phe- 
nomenological triple "ladder" vertex 7. The minus sign in Eq.(4) reflects 
the shadowing origin of this contribution. Eq.(4) can also be derived from 
the calculation of the diffraction dissociation cross section using the AGK 
cutting rules |§, which lead to the relation || 

AF 2 = - F 2 D (5) 

We note that F 2 D , in the above equation, is related to the total integrated 
diffractive DIS cross section. Namely, to the two DIS single diffraction chan- 
nels, as well as the DIS double and central diffraction. It has recently been 
suggested that these are quite large |i0| |TT| . The restriction implied by Eq.(5) 



hinders our ability to reconstruct "theoretically" , the experimental high en- 
ergy behaviour of F 2 (W,Q 2 ). The BFKL approach, which qualitatively re- 
produces the low energy features, fails to do so at high energies. We conclude 
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from this that the necessary corrections, and therefore the diffractive com- 
ponent, must be quite large. As we shall see, when discussing the results of 
our calculation, we require that > 0.25. This requirement is not in con- 
tradiction to the meagre DIS experimental information presently available, 

pD 

allowing one to check Both ZEUS[|r2] and H1|T^] collaborations find a 
sizeable diffractive component in their Q 2 ~ photoproduction studies. In 



DIS, only single diffraction at the 7* vertex has been measured[n][|TjJ, and 
the ratio of the measured diffraction to the total DIS cross section is about 
0.15. Estimates of the non measured diffractive channels are model depen- 
dent. Irrespective of our detailed estimate of the unmeasured channels, the 
overall ratio obtained, is sufficiently large to justify our approach. 



As stated, we take ujq = 0.5. Seemingly, with value of co>o we reach the 
unitarity limit 0], and expect the diffractive channel to vanish. As s-channel 
unitarity is not built into this formalism, we can take u values even larger 
than 0.5, provided that unitarity corrections, such as screening, are incor- 
porated in the calculation. Absorbing all the unknown factors in a new 
phenomenological constant Go we obtain 



W Tr(ln^) 
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F 2 (x B ,Q 2 ) = G ■ {-== ■ exp[ - , OAr J'L w ] (6) 



We have added a (small) constant C to account for the remnant non BFKL 
contributions. 



We now turn to a more detailed discussion of the diagram shown in Fig.2. 
Our motivation is three fold: 

1) We wish to better comprehend the complicated calculation of Bartels, 
Lotter and Wuesthoff and its consequences. 

2) We need to clarify how trustworthy our perturbative QCD calculation is. 

3) We need to adjust the results to the relevant kinematic domain at HERA. 
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We use the expression given for our diagram in Ref . |J 



In— 1 rQ 2 rlk 2 Tr, O 2 
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The integration over k 2 in the above expression leads to an infrared diver- 
gency as k 2 — > 0. However, the BFKL Pomeron is associated with an anoma- 
lous dimension 7(0*) = |, and thus we get Eq.(l) for which the BFKL gluon 
distribution is given by 
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Substituting Eq.(8) in Eq.(7), we rewrite Eq.(7) in a more compact form 
using new variables yp = In—, y B = In—, rn = In and r = In — 

xp xb 

This yields 



af 2 = - 2. r dyp r*^ ( f£,^) . * . e — -«5 (9) 
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where A = 56((3)a s and a s = ^Sa. 

We wish to stress that the infrared divergence of the above integral should 
be studied in more detail. To this end we consider the situation where tq 
and yp are sufficiently large so that we can use the solution of the GLAP 
evolution equation in the region of small to assess F GLAP in Eq.(9). We 
obtain 

F GLAP = A e 2 V Bls( - yB ~ yp ^ r Q^ (10) 

We fix a s so as to perform our calculations in a way consistent with the BFKL 
equation. There is no danger in doing so, as the Q 2 variation in the small 
xp HERA kinematic region is negligible, allowing us to use this approach for 
the analysis of the HERA data. We absorb all irrelevant factors appearing 
before the exponential in a constant factor A, which appears in front of the 
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expression. 



Substituting F,f LAP in Eq.(9) we reduce the equation to the form 
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It is easy to see that there is no saddle point in the integration with respect 
to r. Indeed, the equation for the saddle point is 

where 



i 2r 2 

* = 2Ja a (y B -yp)(r Q -r) + 2uj y P - — — (13) 



Eqs.(12,13) give 
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(14) 



The saddle point can only be at negative values of r, but one cannot trust 
the BFKL equation in this domain, where the virtuality k 2 is less than Q^. 
At such small values of virtualities there are certainly large corrections, and 
it does not seem reasonable to expect the BFKL Pomeron description to be 
valid in this region 0. 

We note that the most important region of integration is still r — > 0, or 
in other words, the dominant value of k 2 remains k 2 ~ Qq. This leads us to 
conclude that the BFKL contribution is questionable and one needs to study 
the integral of Eq.(ll) in more detail, so as to be sure of the domain where 
it is valid. To this end, we observe that our integral over yp has a very good 
saddle point. Indeed, the equation for this saddle point is 
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4 The position of the saddle point in the region of small virtualities has been studied in 
all details in Ref.[p|. 
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Neglecting the last term, we have the saddle point value for yp 



sp a s(rQ - r) 

Vp = vb -r~2 — ( 16 ) 

We now check the reliability of the GLAP approach for the calculation of 
F^ LAP . We recall that the typical value of uj, the argument of the anomalous 
dimension of the GLAP equation^], is given by 



u = J^^A (17) 

V Vb-Vp 

Substituting y P = yp P , we have u = 2uj . 

The BFKL anomalous dimension is given by the series 



7H = ^ + 2C(3) (^) 4 + O(^f) (18) 

UJ UJ UJ° 

Substituting uj = 2uj , we see that the BFKL corrections are very small. This 
does not mean that we do not need the normal GLAP corrections, which are 
essential (see ref. [17|]), but they cannot change the main result of the present 
problem. 

Substituting y P = yp P in Eq.(lO), we end up with the following integral 
over r to be inserted in Eq.(9) 

rQ dre ^V B + ^(r Q -r)-^ (19) 



for ys 3> " s ^3 r ' > - This is the kinematic region which is most interesting 
both from the theoretical and experimental points of view. One can see that 
the integral over r in Eq.(19) is concentrated at small r ~ ^ oc 0(a s ) and 
at large rq and ys- We face the dilemma of how much trust one can put 
on the perturbative calculation which estimates the behaviour of the deep 
inelastic gluon distribution, with virtualities of the order of Q\. Apparently, 
this calculation is not reliable as the problem reduces to that of the energy 

5 We denote u> = N - 1, where N is the moment variable. 
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behaviour of the typical hadron - hadron interaction at high energy, which is 
described by the "soft" Pomeron||T8H [|H| . 



Despite this reservation, the above statement is certainly correct if we 
consider only small a s and/or large Q 2 , the ln(l/xs) parameters in our cal- 
culation. Numerically, the situation is more promising. Indeed, as we shall 
show, the HERA data on F 2 0[[| confirm the theoretical expectation that 
the BFKL Pomeron with u = 0.5 contributes at low xb- Substituting ojq 
= 0.5 we obtain a typical value of r ~ ^ ~ 6, in the integral of Eq.(19). 
This value is sufficiently large to justify our using pQCD, to evaluate the 
diagram of Fig.2. Moreover, r pa 6 is larger than the value of tq in HERA 
kinematic region, so we can estimate the value of the integral in Eq.(19) as 
tq. Collecting all factors together, we obtain Eq.(4), which was used in our 
description of the HERA data. 

Detailed comparisons of our calculations with the data|]]||2| are displayed 
in Fig. 3 where we present Fi(W, Q 2 ) and F2(xb,Q 2 )- We did not attempt a 
"best fit", nevertheless, our ability to reproduce the gross features of the data 
is evident. The following comments relate to the data choice and detailed 
features of our fit: 

1) The data considered are bounded by xb < 10 -2 and W > 50 GeV. 

2) The following parameters were used in the numerical fit: Go = 0.024, 

7 = 0.015 GeV 2 , Q = 1 GeV and C = -0.025. With these parameters we get 
that pa 0.30 at Q 2 = 8.5 GeV 2 and W pa 250 GeV. We can reproduce the 

F D 

data with a which is smaller, but, clearly, our requirement for a relatively 
large DIS diffractive component is essential for this approach. 

3) As can be readily seen from Fig.3, we obtain a reasonable description of 
the data down to values of Q 2 ~ 3 GeV 2 . At smaller values of Q 2 we require 
larger SC to reproduce the data. 

Our inability to reproduce the lowest Q 2 data is not surprising. Clearly, 
the present approach is over simplified, as we fail to take into account the 
different effects of the SC on the diffractive and the total DIS cross sec- 
tions. We call attention to the observation ||1 1|| |T9"|j that in soft hadron in- 
teractions adiff(withSC) pa 0.3(7^// (without SC) whereas a t {with SC) ~ 
0.75a t (without SC). As we have shown [|20f| the strength of the SC in pQCD 
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is determined by a parameter k = ^0z[xG(x,Q 2 )}, where B is the elastic 
slope. At low Q 2 , k ~ 1, and we doubt the validity of Eq.(5), which is at 
the core of our calculation. Nevertheless, even when one can still use 

Eq.(3) with AF 2 as defined in Eq.(4), and obtain a better assessment of F 2 
at low Q 2 than the one we have presented here. The reason for this is, that 
the SC to the total DIS cross section turns out to be smaller than the SC 
to diffraction dissociation. This problem is connected to the broader issue 
of the transition between the (hard) pQCD and the (soft) non perturbative 
domain, which we plan to discuss in a forthcoming publication. 
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Figure Captions 

Fig.l: W dependance of F 2 (W,Q 2 ). Data is taken from ZEUS0 and Hl@ 

with Q 2 > 8.5 GeV 2 . 

v 

Fig. 2: Diffraction dissociation in perturbative QCD. 

Fig. 3: Comparison of F 2 (W,Q 2 ) and F 2 (x,Q 2 ) data with our calculations. 
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